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ABSTRACT 

We present measurements of 68 classical Cepheids with periods from 10 to 78 days 
observed in the near-infrared by the PHAT Program using the Wide Field Camera 3 
(WFC3) on the Hubble Space Telescope (HST). The combination of HST's resolution 
and the use of near-infrared measurements provides a dramatic reduction in the dis- 
persion of the Period-Luminosity relation over the present optical, ground-based data. 
Even using random-phase magnitudes we measure a dispersion of just 0.17 mag, im- 
plying a dispersion of just 0.12 mag for mean magnitudes. The error in the mean for 
this relation is 1% in distance. Combined with similar observations of Cepheids in 
other hosts and independent distance determinations, we measure a distance to M31 
of iiQ = 24.42 lb 0.05{statistical) ± 0.03{systematic), 765 it 28 kpc, in good agreement 
with past measurements though with a better, 3% precision here. The result is also in 
good agreement with independent distance determinations from two detached eclips- 
ing binaries allowing for an independent calibration of the Cepheid luminosities and a 
determination of the Hubble constant. 

Subject headings: cosmology: observations — cosmology: distance scale — galaxies: 
distances and redshifts — stars: variables: Cepheids — cosmological parameters 

1. Introduction 

M31, the nearest analogue of the Milky Way Calaxy, has long provided important clues to 
understand ing the scale of the Universe. T he naked-eye visibility in 1885 of a supernova in M31 



(S And, see lDe Vaucouleurs &: Corwinlll985l ) once suggested the spiral nebulae were located within 



the Milky Way lest its luminosity be "on a scale of magnitude such as the imagination recoils from 
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contemplating^^ (jClerkdll903l ). iHubbld 's (jl929l ) subsequent discovery of Cepheid variables in M31 
revealed the true gulf that existed between the Milky Way and other galaxies. The ability to resolve 
the stellar populations of M31 and knowledge of its distan ce still provides th e best constraints on 
the timescales of the formation of such massive galaxies (IBrown et al.l l2006l ) . Despite long past 
success in identifying Cepheids in M31 (jBaade Sz Swopelll963l . Il965l ). comprehensive inventories of 
its variables were not completed until the availability of wide-format CCD arrays to survey its 3 
degree span and image-subtraction techniques to contend with the extreme crowdi ng of its stars 



a cons equence of its 77° inclination. Large-scale variabilit y surv eys in the 1990 s bvlMagnier et al 



(| 19971 ). the DIRECT Program (e.s. Isonanos et al.l[20oi). and Ivilardell et all (|2007|) as well as 



microlensing surveys in the 2000s, by POINT- AGAPE (I An et al.ll2004l ) and WeCAPP (IFliri et al 
20061 ). succeeded in discovering ~ 10'^ Cepheids, most from > 50 Myr-old supergiants at short 



periods {P < 10 days). The Pixel Observations of M31 with MEgacam (POMME) survey has been 
the most prolific to date, making use of Megacam on CFHT to discover more than 2500 Cepheids 
(jpiiri et ahlboid l. 



The utility of Cepheids as distance indicators critically depends on the accuracy of their mea- 
sured fluxes. Even the best optical measurements of M31 Cepheids from the ground have been 



biase d bright at the 0.1 — 0.2 mag level by stellar crowding (iMocheiska et al. 



200C 



Vilardell et al 



20071 ) and contaminated by variable extinction owing to the high inclination of its host. Space- 
based observations in the near-infrared have the ability to greatly mitigate both of these sources 
of error. A modest sample of 8 Cepheids in M31 was measured with greater resolutio n and in the 



near-i nfrared using the NICMOS Camera on the Hubble Space Telescope (HST) by iMacri et al 



(|200ll ). suggesting that both crowding and extinction could be tamed by such measurements. 



The Panchromatic Hubble Andromeda Treasury (PHAT) Program (P.I. J. Dalcanton) is a 
Hubble Space Telescope Multi-cycle Treasury program to map roughly a third of M31's star forming 
disk, using 6 filters covering from the ultraviolet through the near-infrared. With HST's resolution 
and sensitivity, the disk of M31 will be resolved into more than 100 million stars, enabling a 
wide range of scientific endeavors. Despite capturing only a static view of M31, the random phase 
observations of known Cepheids obtained by PHAT are nearly as precise as mean phase observations 
for meas uring the distance to M31 du e to the small amplitudes of Cepheid light curves in the near- 
infrared (jMadore &: FreedmanI Il99ll ) . Because the Cepheid observ ations are obtain ed with the 
same near- infrared photometric system as recent distance scale data (jRiess et al.ll201ll ). the PHAT 
observations provide the means to determine the distance to M31 with lower systematic error than 
past estimates. The recent discovery and characteriza tion of detached eclipsing binary sytems 
(DEB) in M31 (jRibas et al.l l2005l : IVilardell et al.l l20ld ) offers reliable and independent distance 
determinations to M31 which can be used to calibrate Cepheid luminosities and ultimately the 
Hubble constant. These measurements are also valuable for characterizing the slope of the period- 
luminosity relation at solar me tallicity at long-p eriods in the same near-infrared band used to 



determine the Hubble constant (iRiess et al. 



2011 



We have analyzed the first year of PHAT data to locate and measure the long-period Cepheids 
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(log P (day) > 1) previously discovered by image-subtraction in ground-based variability surveys, 
mostly from POMME. In Section 2 we present images of the recovered long-period Cepheids from 
the PHAT survey and their near-infrared photometry. In Section 3 we analyze this data to constrain 
the distance to M31, the slope of the P~L relation in F160M^ and its impact on the Hubble constant. 



2. WFC3 Observations of Cepheids in M31 



The PHAT program (jPalcanton &: et al.ll2012l ) is imaging the northeast quadrant of M31 with 
WFC3 with 2 UV filters (F275T^, F336W), 2 IR filters {FIWW, F160I^) and 2 optical filters 
with ACS {F475W, i^814iy) over the course of a 3 year survey. All WFC3 images are collected in 
a single epoch with those from ACS obtained approximately 6 months before or after the WFC3 
images when the orientation of HST allows overlapping coverage in parallel. We searched the PHAT 
data for Cepheids with log P > 1 as these match the period range of Cepheids observed by HST 
at > 20 Mpc used in studies of the distance scale. Our primary source of positions and periods 
of Cepheids was the recent POMME Survey from Megacam on CFH T whose sample of > 2500 
Cepheids covering most of M31 is the largest sample collected to date (jFliri et al.ll2012l ). with 180 
in the North half of M31 with log P > 1. While many of the Cepheids found in previous surveys are 
included in the POMME sam ple, the reverse is no t true. We also consulted the variability catalog 
from the DIRECT program (jBonanos et al.l 120031 ) to search for additional long-period Cepheids 
not included in the POMME sample. In the first year PHAT obtained data through the middle 
of 2011, and 65 long-period POMME Cepheids were contained in the survey, a fair fraction of the 
~ 90 likely to be imaged by the end of the Survey. These 65 were a ugmented with 1 object from 
DIRECT and 2 from the Pan-STARRS survey data (|Lee et al.ll201ll ). 



Precise positions of the Cepheids in the HST images of M31 were determined by refining steps 
in the relative astrometry between ground- and space-based imaging. The ground-based survey 
positions were first used to identify the approximate region hosting the Cepheids in the HST 
images to ~ 0'.'4 precision, with the uncertainty resulting from the HST guide star position errors. 
Next, a geometric transformation was defined by matching unresolved sources in 20" diameter 
regions in the ground-based i-band POMME data and the WFC3 FllOVl^ images. The derived 
transformations were used with the POMME Cepheid coordinate to locate the Cepheid positions 
to within one WFC3-IR pixel ((7~0'.'l). For a few Cepheids with close neighbors, identifications of 
the Cepheids in the -F336VF images were used to confirm the Cepheid position among neighboring 
red giants (missing in the UV) or luminous blue dwarves (missing in the IR). Lastly, centroids for 
the Cepheids in the HST images were measured together with photometry resulting in the positions 
given in Table [1] The high signal-to- noise ratio of the Cepheid data in the HST images (SNR > 
100) and minimal crowding (see Figure [T|) insures negligible bias in the measured positions. 
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Table 1. M31 Cepheids in WFC3-IR 
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Table 1 — Continued 
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"Source of optical parameters for the Cepheids: yn*/vs* POMME llFliri et al]|2012h . 0* (Pan-STARRS); D* DIRECT 

llBonanos et a~ ]|20o5^ 
includes absolute zeropoint uncertainty of 0.03 mag 



Table 2. PLCZ Fits to M31 Cepheids 
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The F160W images for each pointing, consisting of 1600 s in 4 dithered exposures, were 
combined and resampled to a final scale of O'.'OS pixel^^. The Cepheid photometry was measured 
by simultaneously fitting model PSFs to the Cepheid and any unresolv ed sources in its vi cinity using 
the same zeropoint scale derived from the standard star P330E as in lRiess et alj (|201ll ). Artificial 



stars were added to the images and measured to assess the crowding bias and to determine the 
photometric errors. The mean bias was 0.002 mag and the mean statistical error was 0.01 mag. 
With the i^llOTy filter, only a single dither of 700 or 800 s was obtained limiting the value of 
resampling the image on a finer scale and the use of PSF fitting, so we measured the FllOVF fiux 
in small apertures of 2 pixel radius. Cepheid parameters are given in Table [TJ 

The F160W and FllOVF P-L relations are shown in Figure [2] with relevant parameters in 
Table [2j For historical interest, we h ave inclu de d one additional measurement in Figure [21 the 
HST WFC3-IR FllOW observation of lHubbld 's k^ld ) first Cepheid discovered in M31, VI with 
P = 31.4 days, observed by his namesake telescope by the Hubble Treasury Program Team (P.I. 
K. Noll). The fitted slope of —3.00 it 0.13 for FIGOVF is in good agreement with the value of 
—2.91 lb 0.06 foun d for 448 more distant Cepheids in 9 hosts with the same filter and instrument 
(jRiess et al.ll201ll ). The M31 slope is about 1.5 a shall ower than the slope of 



3.20 ± 0.06 we 

estimate for the LMC Cepheids from lPersson et al.l (|2004l ) after interpolating between the ground- 
based J- and //-band slopes. The dispersion of the M31 Cepheids about this relation, 0.17 mag, 
is a factor of 3.5 times lower than that of the optical P-L rel ations measured fo r M31 from the 
ground and even less than the 0.23 mag dispersion measured bv lMacri et al.l (j200ll ) with NICMOS, 
likely a result of the greater photometric st ability of WFC3-IR. Although the dispersion is larger 
than the 0.13 mag dispersion measured by iPersson et al.l (|2004l ) in the //-band in the LMC, the 
difference is rea dily explained by ou r use of random phase measurements. Resampling random 
phases from the lPersson et al.l (|2004l ) light curves yields an average dispersion of 0.18 mag with no 
offset, a result nearly identical to ours, and confirming a similar intrinsic dispersion of 0.12 mag. 
A slope-insensitive distance indicator for our sample, the mean Cepheid magnitude at the sample 
mean period of log P = 1.2, is 18.292 it 0.021 mag, sufficient to measure the distance to M31 to 1% 
given sufficient calibration. 

The dispersion increases to 0.20 mag for /^'llOVF. While some of the increase may result from 
additional differential extinction, most appears to come from the larger amplitudes of the light 
curves at shorter wavelengths. The LMC Cepheids predict a random phase dispersion of 0.22 mag. 
A strong correlation between the IR band residuals is apparent in Figure [H as expected from our 
use of random but coincident phases as well as from intrinsic variation. 

Next we fit a Wesenheit relation (PLW in Table [2]) to account for the effect of differential 
extinction along the inclined line of sight of the form F160W — l.BAjFllOW — F160W) =a+61og-P 
where 1.54 is the value of ^FieoVF per magnitude of Apiiow — ^Fwow for a Cardelli et al. (jl989l ) 
reddening law with Ry = 3.1. The slope of this fit, —3.43 it 0.17, is in good agreeme nt with the 
slope of —3.38 lb 0.09 for the same relation in J and H for the LMC (jPersson et al.ll2004l ). A similar 
result was found using the F336W UV color in place of /^'llOVF with a reddening term of 0.14. Next 
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we included a metallicity parameter (PLCZ in Table [2]) to account for any apparent correlation 
between Cepheid near-infrared fluxes and the local value of 12 + \og[0/H] meas ured from HII 



region s, a proxy for Cepheid metallicity. Use of the deprojected radial gradient from lZaritskv et al 



(| 19941 ) results in an insignificant correlation of — 0.65 ±0.73 mag pe r dex, compared to —0.10 it 0.09 
mag per dex from the extragalactic Cepheids in lRiess et al.l (j201ll ). The M31 Cepheids have little 
grasp on the Cep heid metallicity parameter because they lie in a narrow annulus along the disk 
(|Fliri et allboid l with a full range of less than 0.2 dex (12 + log[0/i7]=8.87 to 9.05). 



To determine the distance to M31 as well as other parameters of int erest we now make use of 
the sample of Cepheids measured in the near-infrared bv lRiess et al.l (j201ll ). multiple anchors for the 
Cepheid distance scale, an d SN data which ca n be used to determine the Hubble constant. We follow 

solving one simultaneous system of linear equations 



the same formalism as in iRiess et al.l (2011 



which relate Cepheid magnitudes to their absolute magnitudes, slope of their P-L relations, local 
metallicity, distances to their hosts, and -with the inclusion of SN la data- a measurement of the 
Hubble constant. The Cepheids of M31 are assumed to share the same nuisance parameters as other 
Cepheids (i.e., luminosity, slope of P-L , metallicity dependence) but with a unique distance. For 
the M31 Cepheids, the photometric system used to measure their colors was somewhat different. 
While the Cepheids in the 8 SN la hosts, the maser host NGC 4258, and M31 were all measured 
with F160W on WFC3-IR, the optical colors of the POMME M31 Cepheids, useful for dereddening, 
were not measured by Megacam or PHAT with the same V and / bands on HST as the others. To 
account for this difference we employed one of two different prescriptions: 1) we assumed a uniform 
value for the V — I color of the M31 Cepheids as the mean of those measured from the ground by 
the DIRECT program, V-I = 1.23±0.03, (indicated as fit PLW = Hyj in TableE]); or 2) we used 
the individual FllOT^ — -F160VF colors measured for the M31 Cepheids from the PHAT data with 
a small offset derived to give the same mean color correction in V-I from DIREClQ (indicated 
as PLW = Hj^H,x in Table [3]). The advantage of the latter approach is that it can account 
for differential reddening along the line of sight while providing a reddening correction which is 
consistent with that used for non-M31 Cepheids. We adopt an 0.03 mag systematic uncertainty 
for the use of colors measured with a different photometric system and an 0.04 mag systematic 
uncertainty between near-IR magnitudes of Cepheids measured on the ground and those measured 
from space. 



^By equating the mean V — I dereddening with that for FllOW — F160W we can solve for a color offset to insure 
they have the same mean. That is, 0.504 < V - I >= 1.54 < FllOW - F160W -X >, where <V-I >=1.23, gives 
X = -0.066 mag. 
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° Note: metallicity calibration reference: zkh: IZaritskv et al.l (|19941 ). Tg: iBresolinI (|201l[ ). 
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Our best distance estimate for M31 is fiQ = 24.415 it 0.052, a 2.4% (statistical) distance 
determination which makes use of independent distance determinations t o NGC 4258, Milky Way 
Cepheid parallaxes, and DEB distances in the LMC (see lRiess et al.ll201ll . for a description of these 
distance-scale anchors) . 

The best fit parameters are given in Table EJ top line. Column 2 gives the value of Xdof^ 
Column 3 the distance modulus of M31 and its statistical uncertainty, Column 4 the number of 
Cepheids used in the fit, Column 5 the value and total uncertainty in Hq, Column 6 is a flag to 
indicate the use of Cepheids below the optical completeness limit (see below). Column 7 gives the 
metallicity calibration used. Column 8 the correlation coefficient in the PLWZ regression. Column 
9 the value and uncertainty of the slope of the Cepheid PLWZ relation. The next three columns 
are used to indicate variants in the analysis whose impact we now consider. 

Ignoring the metallicity parameter (lines 2 and 4), changing the method of color correction 
(line 3), changing the reddening paramete r for the Ceph e ids fr om Ry = 3 .1 to Ry = 2.5 (line 



6), or changing the metallicity scale from IZaritskv et al.l (| 19941 ) to that of iBresolinI (|201ll ) (line 
8), changes the distance to M31 by < 0.02 mag. Failing to clip outliers in any Cepheid P-L 
relations (line 7) doubles the Xdof changes the distance to M31 by < 0.01 mag. Retaining 
infrared measurements of Cepheids with periods below the optically-determined completeness limit 
(indicated by < P = A^, line 8), has no effect on the distance. The only significant change in the 
distance to M31 occurs when discarding the color measurements used to account for extinction 
(line 5). The 0.15 mag increase in the distance in this case indicates that the M31 Cepheids have 
more extinction, AAh = 0.15 mag (or DeltaAy = 0.75 mag) than the average of the mostly 
extragalactic Cepheids in other hosts. This is not surprising as M31 is more inclined than any of 
the other hosts. We also used each of the 3 distance-scale anchors separately as shown in lines 10 
to 12. The uncertainties in the distances increase with the us e of only a sing l e anch or. The lowest 
uncertainty is from the use of the parallax measurements by [Benedict et al.l (|2007l ) to Milky Way 
Cepheids (line 11), resulting in ixq = 24.38 ± 0.064. 



Following the approach of iRiess et al.l (|201ll ). we quantify the systematic uncertainty in the 
distance to M31 from the dispersion of the variants in the analysis, a = 0.03 mag for the 12 variants 
not including the one neglecting color information which ignores the large extinction for the M31 
Cepheids due to inclination. Thus our best estimate of the distance to M31 is fiQ = 24.42ib0.05±0.03 
or 76 5 ifc 28 kpc. This dis t ance is in excellent agreement with the frequently-cited measurement 
from iFreedman &: Madord (Il99d) and sit s near the middle of the range of past measurements 
summarized bv iMcConnachie et al.l (|2005l ). 



Lastly, we made use of the two DEB measurements for M31 from Ribas et al. (j2005l ) and 
Vilardell et al.l (j2010l ) with a mean of = 24.36 it 0.08 to determine the Hubble constant by their 
ability to calibrate the Cepheid luminosities. Combined with the preceding three anchors (line 13), 
the use of the independent M31 distance has negligible impact, increasing the Hubble constant by 
0.2 km s^^ Mpc^^ to 74.0 km s~^ Mpc~^ with no reduction in uncertainty. The present limitations 
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of the use of M31 as an independent anchor are the lack oiV — I colors measured with HST WFC3 
and the lower precision of its independent distance compared to the other 3 anchors. The use of 
M31 DEB results without any of the other anchors (line 14) yields a larger Hubble constant with 
larger error though still consistent with prior results. The best fit to the Hubble constant (line 1), 
which only makes use of the M31 Cepheids to constrain the slope and metallicity p arameter, is 



73.6ifc 2.4 km s ^ Mpc ^ and does not bring a significant difference to the one inferred by lRiess et al 
(1201 ih. 



3. Discussion 

The main advantages of the PHAT space-based observations of Cepheids presented here over 
previous data come from the reduction in extinction and crowding. The result is the tightest 
P-L relation for M31 Cepheids yet seen, which, combined with past, external calibrations, yields 
the most precise distance measurement for M31. The improved precision over the past history of 
optical, ground-based P-L relations is striking and contin ues to indicate th e value of this kind of 



data for measuring the Hubble constant and dark energy (iRiess et al.ll201ll ). 



To see more clearly the advantages of these spaced-based observations over those from the 
ground we simulated the effect of crowding for near-infrared Cepheid magnitudes obtained with 
good ground-based seeing by measuring the flux contained in an aperture of radius r=Cf.'9. The 
dispersion increased from 0.20 to 0.24 mag in i^llOVF, the slope of the P-L relation flattened by 0.5 
and the intercept became brighter by 0.3 mag. This is not unexpected as the impact of crowding 
is greater for the lower-period, fainter Cep heids. These results ar e similar to those determined 



for ground-based crowding in the optical by iMochejska et al.l (120001 ) . It is important to note that 



the effect of crowding is far more severe for M31 due to its large inclination than for less inclined 
galaxies, even those much farther away. In addition, Cepheids found from the ground via image 
subtraction are likely to suffer significantly greater crowding than those selected from PSF fitting 
as the addition of comparable constant fiux to a PSF fit will reduce the amplitude of the light curve 



and remove the object from an amplitude-selected sample (jFerrarese et al.ll2000l ). 



The P-L relations for M31 may still improve in the near future. Additional observations by the 
PHAT program should augment the Cepheid sample by dozens. Cepheids with logP < 1, though 
less useful for distance scale work, can be mined from the data to study the short period end of the 
relation. The use of Cepheid phase information from concurrent, ground-based optical monitoring 
of M31 can be used to recover the phase of the PHAT observations, reducing the scatter of the P-L 
by up to 50%, an equivalent leverage as doubling the sample of measurements presented here. 

We are grateful to the members of the PHAT MCT Program led by Julianne Dalcanton 
and aided by Jason Kalirai for their tremendous efforts to obtain the PHAT measurements and 
their support of this work. Financial support for this work was provided in part by the POM- 
MME project (ANR 09-BLAN-0228). This work is also based on observations obtained with 
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MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the CFHT which is op- 
erated by the National Research Council (NRC) of Canada, the Institut National des Sciences de 
I'Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University 
of Hawaii. JF acknowledges the hospitality of GEPI, where part of the POMME work was carried 
out. 
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Fig. 1. — HST near-IR FIGOVF images of 68 Cepheids in M31 (ordered as in Table [T]), spanning 
nearly an order of magnitude in period. The scale of each stamp is 2.5" (32 pixels). The position of 
the Cepheid as determined from the optical Megacam POMME images in indicated by the circle, 
which has a diameter of 1". 
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Fig. 2. — Near-IR P-L relations for 69 Cepheids in M31 with logP > 1 (Table [T]). The single slope 
fitted to the relations is given in Table [21 and is shown as the solid lines. Dashed lines indicate the 
average dispersion of 0.17 mag (F160VF), a factor 3.5 smaller than previous ground-based optical 
P~L relations, and 0.20 mag (FllOVF). 



